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THEORY  OF  CROSS  FIELD  DEVICES  AND  A  COMPARATIVE 
STUDY  OF  OTHER  RADIATION  SOURCES* 

I.  Introduction 

Magnetron  is  a  widely  used  radiation  source  [e.g.,  radar,  microwave 
oven,  etc.,].  It  was  developed  very  early  but  is  perhaps  the  most 


complicated  one  to  analyze  theoretically 


Design  studies  have  largely 


relied  on  empirical  experiences  and  on  the  two  formulas,  the  Hull  condition 
and  the  Buneman-Hartree  condition,  which  were  derived  almost  half  a  century 
ago.  Modern  advances  in  electronic  computing  still  have  not  provided 
sufficiently  accurate  models  to  confidently  predict  or  explain  the 

multiple  facets  usually  observed  in  experiments.  Such  a  state  of  imperfect 
understanding  of  the  magnetrons  remains  in  spite  of  the  strong  efforts  of 
the  many  leaders  of  the  field.  Resurgence  of  interest  in  cross  field 
devices  is,  in  part,  triggered  by  the  more  recent  developments  of  high 
power  diodes^  and  relativistic  magnetrons^’ ^ . 

Much  of  the  difficulty  in  a  theoretical  study  of  magnetrons  lies  with 
.  their  complex  geometries.  Controversy  already  emerged  at  the  very  first 
step:  What  should  be  taken  as  the  basic  motion  of  the  electrons  in  a 
smooth  bore  magnetron,  in  the  absence  of  rf  fields?  Should  one  take  the 
highly  idealized  model  of  laminar  flow  or  should  one  really  need  to  use  the 
more  complicated  "multi-stream"  model,  in  which  the  electrons  follow 
cycloidal  orbits  which  commence  and  finish  on  the  cathode  surface?  To  this 
date,  a  definite  answer  still  has  not  emerged,  but  general  opinion  seems  to 
be  more  in  favor  of  the  laminar  flow  model,  especially  prior  to  full 
oscillation.  Note  that  the  underlying  words  above  already  indicate 
considerable  simplification.  In  fact,  in  virtually  all  existing  analytical 

*  I  his  m.imisvnpi  is  ,  p  if  c>l  tuf  <>f  Hivh  /*•'»»»•»  \/ ,  »nm  in  xtml  M'llmuh'  Hwm  l  vl\  V  I  in  111  ihK'in  m»|  I  \  \  1 
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theories,  the  presence  of  corrugated  structure  is  assumed  to  have  no  effect 
on  the  unperturbed  motion  of  the  electrons^.  The  wall  corrugation  enters 
only  in  the  rf  modes  which  are  supported  [Fig.  la).  It  is  possible  that 

such  an  assumption,  made  for  convenience  of  analysis,  has  already 
eliminated  a  significant  amount  of  interesting  physics.  [One  needs  not  go 
very  far  to  see  this.  For  example,  it  is  the  periodicity  introduced  in  the 
streaming  motion  of  the  electrons  in  the  equilibrium  state  (by  a  magnetic 
wiggler)  which  is  responsible  for  free  electron  laser  action  as  is 


explained  elsewhere  in  this  volume.) 

Another  feature  which  contributes  to  the  complexity  of  magnetrons  is 
the  fact  that  the  electron-emitting  cathode  is  part  of  the  interaction 

g 

circuit.  This  differs  from  many  other  microwave  devices  such  as  gyrotron  , 
klystron^,  free  electron  laser^,  gyromagnetron^,  peniotron**,  etc.,  where 
the  interaction  region  is  separated  from  that  of  beam  generation.  The 
control  of  gun  noises  and  "beam  quality"  becomes  problematical.  This,  of 
course,  is  just  the  reason  why  the  beam  cannot  be  well-characterized  as 


described  in  the  preceding  paragraphs.  This  problem  of  equilibrium 


modelling  appears  to  be  an  intrinsic  difficulty  for  theorists,  in  those 

microwave  sources  where  the  rf  draws  directly  from  the  electrostatic  energy 

of  the  system.  |Another  case  in  point  is  the  more  recent  invention  of 
12 

orbitron  ,  where  the  unperturbed  orbits  cannot  be  controlled,  thereby 
leading  to  a  variety  of  intei pretat ion  and  interaction  theories.  See  below 
for  more  detail.)  Actually,  this  may  well  be  ar:  advantage  as  far  as  the 
utility  of  these  tubes  is  concerned.  Such  a  tube  should  be  rugged  when 
operational,  since  the  peifoimance  is  not  ciitically  dependent  on  whether 
the  beam  has  a  superb  beam  c|irality.  This  insensitivity  to  the  electron 
beam  quality  is  i n  marked  contrast  to  the  more  "delicate"  radiation  smncc. 
such  as  travelling  wave  tubes,  free  electron  lasers,  etc. 


* 


The  magnetron  circuit  supports  a  variety  of  electromagnetic  modes. 
The  multiplicity  of  modes,  together  with  the  uncertainties  in  the 
electronic  motion,  greatly  reduces  the  predictability  of  beam-circuit 
interaction.  In  fact,  whether  the  beam-circuit  interaction  analysis,  using 
the  standard  technique  of  impedance  matching,  is  adequate  at  all  in 
describing  magnetron  operation,  has  been  questioned  from  the  beginning. 
This  concern  was  not  without  basis  and  was  reviewed  by  J.  M.  Osepchuk  on  p. 
275  of  Ref.  2.  It  perhaps  explains  why  the  theory  of  magnetron  lags  behind 
those  of  the  more  modern  devices  mentioned  earlier.  But  it  also  means  that 
the  full  potential  of  high  power  radiation  sources  driven  by  direct 
conversion  (from  potential  energy)  remains  to  be  realized. 

With  this  caveat,  we  shall  only  give  here  an  elementary  exposition  of 
some  aspects  of  cross  field  devices.  The  laminar  flow  model  will  be 
adopted.  The  operating  conditions  will  be  given  in  the  next  section,  to  be 
followed  by  a  qualitative  description  of  the  mechanism  of  energy  transfer 
in  cross-field  geometries.  The  r£  generation  is  a  result  of  the  interplay 
of  wave-particle  synchronism,  phase  focusing,  conversion  of  both 
electrostatic  and  kinetic  energy  of  the  beam  through  its  transverse  and 
longitudinal  motion  in  the  rf  fields  supported  by  the  corrugated  walls. 
Some  (but  not  all)  of  these  features  are  shared  by  other  more  contemporary 
microwave  sources  such  as  orbitron,  gyrotron,  peniotron  and  gyro-magnetron 
and  will  be  discussed  in  Sec.  III.  In  fact,  the  comparison  among  these 
various  devices  form  an  integral  part  of  this  Chapter.  This  can  be 
achieved  by  the  use  of  highly  idealized  models  which  will  be  described  in 
Sec.  IV  and  V.  Despite  its  scope,  the  exposition  is  self-contained,  and  a 
minimal  amount  of  mathematics  will  be  required.  Detailed  justification  and 
extension  of  the  theory  will  be  referred  to  when  the  occasion  arises.  Some 
of  the  topics  discussed  are  still  under  active  investigation. 


II.  Brillouin  Flows,  Hull  Cut-off  Condition,  and  Buneman-Hartree 

Oscillation  Condition 

The  laminar  flow  model  which  has  been  studied  extensively  in  cross 
field  devices  is  the  Brillouin  flow.  It  is  basically  an  ExB  drift, 
ignoring  the  cycloidal  motion  of  the  electrons  [Fig.  1].  It  is  the  mean 
flow  field  of  the  electron  cloud  with  which  the  rf  interacts  over  a  long 
time  scale.  Although  the  Brillouin  flow  is  difficult  to  achieve 
experimentally,  recent  particle  simulation  of  magnetrons  show  that  the 
space  charge  distribution  and  the  electronic  motion  contain  a  large 
component  of  the  Brillouin  flow,  especially  during  the  start-up  stage.  In 
the  absence  of  curvature  effects,  consideration  of  the  planar  Brillouin 
flow  immediately  leads  to  the  Hull  cut-off  condition  and  the  Buneman- 
Hartree  condition.  These  conditions  determine  the  range  of  voltage  (at  a 
given  magnetic  field)  in  which  a  magnetron  operates.  When  curvature  and 
relativistic  effects  are  included,  the  derivation  would  be  more  involved, 
but  the  main  features  of  these  conditions  remain  unchanged. 

One  should  keep  in  mind  the  assumptions  used  in  the  derivation  of  the 
Brillouin  flow:  (a)  The  flow  is  laminar,  and  is  parallel  to  the  sole. 

(b)  There  is  no  variation  in  the  direction  of  the  flow.  Thus,  spatial 
inhomogeneity  associated  with  the  corrugated  structure  is  ignored. 

(c)  All  electrons  are  emitted  from  the  cathode  with  zero  velocity  and  zero 
energy.  In  the  derivation  of  the  cut-off  condition  and  the  oscillation 
condition,  the  rf  fields,  when  present,  are  assumed  to  have  infinitesimally 
small  amplitude,  and  thus  all  fields  are  essentially  static.  (The  last  two 
assumptions  are  also  used  when  these  conditions  are  derived  from  single 
particle  motion  as  in  the  case  of  a  cylindrical  model.)  The  generalization 
to  cylindrical  models  with  relativistic  effects  will  be  indicated. 


For  simplicity,  we  shall  first  consider  the  planar,  nonrelativistic 


model.  The  governing  equations  for  the  system  depicted  in  Fig.  lb  read 


v  =  x  v(y) 
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Equation  (1)  follows  from  assumptions  (a),  (b),  and  (2)  from  assumption  (c) 
of  the  previous  paragraph.  Equation  (3)  accounts  for  the  force  balance  in 
the  flow  field.  Here  v(y)  is  the  flow  speed  and  V=V(y)  is  the  voltage 
distribution  with  V(0)=0  at  the  cathode,  e  is  the  electronic  charge,  Eq  is 
the  electric  field  consistent  with  the  space  charges  and  Bq  is  the 
external  constant  magnetic  field.  The  diamagnetic  field  is  ignored. 

Upon  differentiating  Eq.  (2)  with  respect  to  y  and  eliminating  3V/3y 
according  to  (3),  we  obtain  a  solution  for  v(y): 

v(y)  =  <»cy  (4) 

where  «c  =  |e|Bo/mo  is  the  cyclotron  frequency.  In  obtaining  (4),  we  have 
used  the  boundary  condition  v(0)=0.  The  voltage  distribution  reads 


V(y) 
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and  the  electric  field  is 
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The  charge  density  is 
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which  may  be  written  as 
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where  »  =  (p  e/m  e  )  is  the  the  plasma  frequency.  Note  that  the 

p  o  o  O  j 

laminar  flow  condition  demands  the  condition  w  =co  be  satisfied. 

p  c 

Given  separation  D,  operating  voltage  V  and  magnetic  field  Bq,  the 

thickness  H  of  the  space  charge  layer  {Fig.  lb]  may  be  calculated.  Note 

2 

first  that  the  electric  field  equals  to  -mo<*)^H/|e|  between  the  region  y=H 
and  y=D  [cf.  Eq.  (5)].  Then  the  total  voltage  drop  becomes 


2U 
m  w  H 
o  c 


V  =  V(H)  +  <D  -  H>  =  “cBoH(D  "  V' 


(7) 


Solving  for  H,  we  obtain 

H  =  D [l  -  (1  -  V/Vc)1/2],  (8) 

where 

V  =  \  B2D2  (9) 

c  2  m  o  v  ' 

o 

is  the  Hull  cut-off  voltage. 

Thus,  the  Brillouin  flow  is  characterized  by  a  velocity  field  of 

constant  shear  [cf.  Eq.  (4)],  and  by  a  constant  charge  density  extending 

from  the  cathode  to  a  distance  H,  which  is  given  by  (8)  in  terms  of  the 

voltage  V,  magnetic  field  Bq  and  anode-cathode  separation  D.  The  electron 

density  satisfies  w  =w  . 

pc 

The  Hull  cut-off  voltage  V  has  the  following  interpretation.  It  is 
the  maximum  voltage,  at  a  given  magnetic  field,  allowed  for  magnetic 
insulation,  i.e.,  for  prohibiting  the  electrons  from  reaching  the  anode. 
This  can  be  seen  by  setting  H  =  D  in  Eq.  (7).  The  parabolic  dependence  of 
Vc  on  Bq  [cf.  Eq.  (9)]  is  sketched  in  Fig.  2.  For  parameters  below  t  e 
Hull  parabola,  the  magnetron  is  insulated.  The  same  condition  (9)  is 


obtained  from  single  particle  orbit  considerations.  An  electron  emitted 

from  rest  at  the  cathode  will  not  be  able  to  reach  the  anode  if  V<V  ,  if 

c 

the  fields  are  time  independent  and  azimuthally  symmetric. 

An  alternate  form  for  Eq.  (8)  is 


(  r 

oqi/2) 

H  =  D  I1  -  [ 

1  -  <VBo>  ]  j 

where 


BH 


1/2 


(10) 


is  known  as  the  Hull  cut-off  magnetic  field.  It  is  the  minimum  magnetic 
field  required  for  magnetic  insulation  at  a  given  voltage.  When  the 
operating  voltage  V  is  small  compared  with  V  ,  we  may  approximate  (8)  as 


which  shows  that  the  space  charge  layer  thickness  increases  linearly  with 

2 

the  voltage,  but  decreases  as  1/Bq  as  Bq  increases.  If  the  operating 
voltage  is  only  1/2  of  the  Hull  cut-off  voltage,  the  space  charge  layer 
occupies  about  1/4  of  the  width  of  the  gap.  It  is  easily  shown  from 
Eq.  (8)  that  as  V  -»  ,  magnetic  insulation  is  rapidly  diminished. 

For  effective  transfer  of  energy  from  the  electrons  to  the  rf,  there 
must  be  synchronism  between  the  rf  wave  and  at  least  some  of  the  electrons 
in  the  sheath.  Since  the  ExB  drift  of  the  electrons  within  the  sheath  has 


a  maximum  speed,  for  rf  with  phase  velocity  v  ^  to  be  synchronized  with  the 


electron  at  a  fixed  magnetic  field,  the  voltage  cannot  be  too  low.  For 
then  the  space  charge  layer  thickness  H  would  be  small,  and  even  the 


maximum  electron  speed  o^H  would  be  below  v  ^  [cf.  Eqs.  (4),  (8)].  This 


minimum  voltage,  above  which  a  magnetron  may  oscillate,  is  called  the 
Buneman-Har tree  threshold  voltage  and  is  determined  in  the  present  case 
simply  by 


7 


.TTjT 


v<H)  =  vph.  (11) 

Upon  using  (8)  and  (4)  in  (11),  the  threshold  voltage  then  becomes 

m  _  , 

V_„  =  B  Dv  ,  -  sA  v  f,  (12) 

BH  o  ph  2  { e  (  ph’  v  ' 

which  depends  on  the  phase  velocity  of  the  mode  under  consideration.  This 

threshold  voltage  is  also  shown  in  Fig.  2.  It  can  easily  be  shown  that, 

regardless  of  the  phase  velocity,  the  straight  line  representing  the 

threshold  voltage  in  Fig.  2  is  always  tangent  to  the  Hull  parabola. 

When  both  relativistic  and  cylindrical  effects  are  included,  the  Hull 

cut-off  magnetic  field  and  the  threshold  voltage  are  modified.  If  we  use 

2  2 

D  to  represent  the  equivalent  gap  width  (b  -a  )/2a,  [cf.  Fig.  la],  the 


Hull  cut-off  magnetic  field  becomes 


2,1/2 


H  -  1 


and  the  threshold  voltage  (12)  is  modified  to  read 


(10a) 


me"  (  ,  .v,.2,l/2 

FT  j1  “  t1  -  H11)  1 
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(12a) 


This  relativistically  correct  threshold  voltage  was  obtained  by  Walker 

using  a  Hamiltonian  description  of  orbits.  When  this  threshold  voltage  is 

maximized  with  respect  to  the  phase  speed,  the  cut-off  voltage  [Eq.  (10a)] 

is  obtained.  Again,  the  threshold  voltage  curve  is  tangent  to  the  cut-off 

voltage  curve  as  in  Fig.  2.  Magnetron  oscillation  takes  place  when  V  and 

Bq  lie  between  these  curves,  but  usually  closer  to  the  Buneman-Har tree 

threshold  voltage^.  Note  that  (10)  and  (12)  are  recovered  from  (10a)  and 

(12a)  in  the  nonrelativistic  limit  c  -»  ro.  There  are  recent  modifications 

14 

of  these  voltages,  including  the  effect  of  diamagnetic  current  and  axial 


8 


current*^.  All  of  these  derivations  made  three  key  assumptions:  (a)  the 
external  field  is  azimuthally  symmetric,  (b)  the  amplitude  of  rf  is  so  low 
that  the  energetics  of  orbital  motion  is  not  affected,  and  (3)  the  fields 
are  time  independent.  In  practice,  none  of  these  assumptions  are  valid. 
We  refer  to  the  Chapter  by  Benford*^  for  a  detailed  discussion  of  the 
experiments  and  these  conditions. 

The  validity  of  the  old  concepts  outlined  above  has  been  tested  in 

computer  simulations.  A.  Palevsky^  recently  developed  a  two-dimensional, 

time-dependent,  fully  electromagnetic  and  fully  relativistic  particle 

simulation  code  to  model  pulsed  relativistic  magnetrons*^.  Some  of  the 

simulation  results  are  shown  in  Figs.  3  and  4.  Figure  3  shows  the  contour 

levels  of  charge  densities  at  an  early  time,  before  oscillation  starts. 

The  arrows  label  the  theoretical  Brillouin  layer.  The  corresponding 

velocity  profile  (not  shown)  follows  the  Brillouin  flow  to  within  a  few 

percent.  The  contour  plot  at  a  much  later  time  (when  oscillation  is  fully 

developed)  is  shown  in  Fig.  4.  The  spoke  formation  is  evident.  Note  that 

the  spatial  extent  of  the  sheath  is  approximately  given  by  the  Brillouin 

layer  even  at  this  late  stage  when  the  beam-circuit  interaction  is  highly 

nonlinear.  The  velocity  profile  corresponding  to  Fig.  4  contains 

considerable  spread  about  the  theoretical  Brillouin  flow,  and  confirms  the 

18 

considerable  back-bombardment  of  electrons  at  the  cathode  .  These 

simulation  results  do  show  that  the  electron  flow  is  a  mixture  of  Brillouin 

flow  and  a  somewhat  weaker  multi-stream  component.  They  also  corroborated 

19 

the  suggestion,  made  by  Buneman  a  long  time  ago,  that  the  instabili ties 
associated  with  such  shear  flows  may  play  a  key  role  in  the  start-up  of 
magnetron  oscillation.  We  shall  return  to  the  stability  of  the  laminar 
flow  later. 
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Ill .  Gain  Mechanism  in  Cross  Field  Devices  and  Comparison  With  Other 

Radiation  Sources 

In  this  section,  we  describe  in  a  qualitative  manner  the  transfer  of 
electron  energy  (kinetic  and  potential)  to  the  rf.  We  shall  bring  out 
certain  features  which  are  unique  to  the  magnetron  configuration.  We 
shall  also  contrast  several  other  microwave  sources  which  share  some,  but 
not  all,  of  these  features.  A  simple  theory  using  a  beam  model  which  is 
often  adopted  for  all  of  these  radiation  sources  will  be  given  in  the  next 
section. 

If  the  externally  imposed  voltage  and  axial  magnetic  fields  are 

sufficiently  strong,  the  average  motion  of  the  electrons  is  approximately 

the  ExB  drift.  The  curvature  effect  of  the  magnetrons  is  secondary  in 

importance,  and  one  may  consider  a  planar  model  as  shown  in  Fig.  5.  Let  us 

consider  the  action  on  electron  A  due  to  the  combined  DC  fields  E  ,  B  ,  and 

o  o 

an  rf  electron  field  E^.  The  rf  electric  field  is  assumed  to  have  a  phase 
velocity  approximately  equal  to  the  unperturbed  drift  velocity  vq  of  the 
electron,  and  to  have  a  phase  shift  of  rt  from  one  cavity  to  the  next,  as 
shown  in  Fig.  5.  The  rf  electric  field  Ej  affects  the  motion  of  electron  A 
in  four  ways  :  (a)  It  retards  the  drift  motion  of  the  electron  because  E^ 
and  vq  are  in  the  same  direction.  (A  carries  a  negative  charge).  (b)  The 
EjXBq  drift  brings  test  electron  A  upward,  to  a  location  closer  to  the  wall 
corrugation  (anode),  where  the  rf  field  is  stronger.  (c)  As  a  result,  the 
drift  motion  is  retarded  further,  and  by  conservation  of  energy,  this 
energy  is  converted  to  rf.  (d)  More  importantly,  as  electron  A  drifts 
toward  the  anode,  it  experiences  a  loss  of  potential  energy  associated  with 
the  DC  electric  field.  Again,  conservation  of  energy  requires  that  the 
potential  energy  be  converted  to  the  rf.  These  processes  continue  until  A 


hits  the  anode. 


Consider  now  the  action  of  the  rf  electric  field  on  test  electron  B 
[Fig.  5J.  Test  electron  B  is  accelerated  by  the  field.  But  this 
acceleration  is  weakened  gradually  since  the  E^xBo  drift  brings  electron  B 
downward,  toward  the  cathode  where  the  rf  electric  field  is  weaker.  Thus, 
while  electron  B  gains  energy  from  the  rf,  on  the  whole,  the  amount  of 
energy  it  gains  is  less  than  the  rmount  of  energy  electron  A  would  lose  to 
the  waves,  and  the  rf  signal  grows.  Note  that  this  rf  growth  is  at  the 
expense  of  both  the  kinetic  and  potential  energy  of  the  electron,  the 
potential  energy  giving  the  dominant  contribution. 

The  configuration  shown  in  Fig.  5  also  demonstrates  the  interesting 
property  of  phase  focusing  in  magnetrons.  Suppose  now  the  phase  speed  v  ^ 
of  the  rf  is  slightly  less  than  the  drift  speed  vq  of  the  electrons.  At  a 
later  time,  electron  A  will  advance  ahead  of  the  wave,  reaching  a  position 
near  P  relative  to  the  wave.  Near  P,  opposes  Eq.  Therefore,  the 
electron  drift  velocity  (Eq+E^)xBo  is  lowered,  thereby  reducing  the 
difference  between  v  ^  and  the  drift  velocity.  This  effectively  brings  the 
electron  from  position  P  back  to  the  original  position  A.  For  electron  B, 
since  vq  >  v^,  at  a  later  time,  it  reaches  position  near  Q  with  respect  to 
the  wave.  Near  Q,  E^  is  in  the  same  direction  of  Eq,  therefore,  the 
combined  drift  velocity  (E^+Eq)  x  Bq  is  increased,  enlarging  the  difference 
between  the  drift  velocity  and  the  phase  speed  of  the  wave.  In  othei 
words,  electrons  like  B,  which  gain  energy  from  the  rf,  are  quickly  put  out 
of  synchronism,  whereas  electrons  like  A,  which  give  energy  to  the  rf,  will 
continue  to  do  so  until  they  reach  the  anode. 

The  same  argument  applies  when  v0^vph"  For  example,  if  <  Vph’ 


electron  A  will  lag  behind  the  wave,  moving  *o  a  position  near  ()  witli 
respect  to  the  rf.  At  0,  the  drift  speed  (E.+E  )xB  increases  since  E. 


reinforces  Eq  there.  This  effectively  brings  A  closer  to  synchronise  with 
the  wave  than  was  originally. 

From  the  description  of  the  last  two  paragraphs,  one  sees  that  there 
is  a  built-in  phase  focusing  for  the  "favorable"  electrons,  (like  A  which 
gives  its  energy  to  the  rf)  and  a  phase  defocussing  of  the  "unfavorable" 
electrons,  (like  B  which  gains  energy  from  the  rf).  It  is  this  built-in 
focusing  property  which  is  largely  responsible  for  the  high  efficiency  of 
magnetrons.  It  is  the  reason  why  charge  spokes  are  formed  |cf.  Fig.  4|  and 
is  perhaps  a  physical  basis  for  the  construction  of  nonlinear  soli  ton 
solutions'^  which  are  shape-preserving. 

The  rather  fascinating  interplay  between  DC  and  AC  electric  fields, 
drift  motion  and  the  presence  of  corrugated  structure  in  a  cross  field 
configuration  can  be  contrasted  with  other  microwave  sources  which  shaie 
some  of  these  features.  Ve  have  in  mind  the  peniotron,  orbitton,  gyiotron 
and  gyromagnetron.  All  of  these  devices  are  currently  under  active  study. 
They  were  invented  more  than  twenty  years  after  the  magnetrons  were  put  in 
service.  In  a  sense,  they  are  much  "simpler"  devices. 

A  model  of  peniotron1*  is  shown  in  Fig.  6a.  It  consists  of  an  annular 
electron  layer  rotating  inside  a  ridged  waveguide-  This  beam  is  guided 
only  by  an  external  magnetic  field.  This  system  differs  from  a  magnetron 
in  that  the  rf  energy  derives  from  the  kinetic  energy  of  the  electron  beam, 
the  potential  energy  is  entirely  absent.  (No  DC  electric  field  is  present 
in  the  interaction  circuit.)  It  is  similar  to  the  magnetron  in  otto 
important  aspect,  however.  The  spatial  i nhomogene i t y  of  the  if  fields 
produced  bv  the  ridged  waveguides  may  cause  the  guiding  center  of  an 
election  to  drift  sideways,  and.  as  a  result,  to  turn  virtually  all  of  its 
kinetic  energy  into  radiation,  implying  an  extremely  high  efficiency. 


To  see  the  action  of  the  rf  field  inhoeogenei ty  on  the  energy 

transfer,  consider  two  electrons  A  and  B  in  Figs.  6b  and  6c.  Assume  that 

the  electrons  are  "nonrelat i vist ic" ,  (~1()  KeV,  as  typical  in  peniotron)  and 

that  its  cyclotron  frequency  is  1/?  of  the  rf  frequency  (i.e.,  operation  at 

second  cyclotron  harmonics).  Electron  A  is  in  an  accelerating  phase,  its 

rotational  energy  increases  and  as  a  result,  its  radius  enlarges.  As  its 

radius  enlarges,  it  aoves  half  a  cyclotron  period  later  to  a  position  whete 

the  rf  field  is  stronger,  but  at  which  time  the  if  field  decelerates 

electron  A.  Thus  electron  A  gives  moie  energy  to  the  if  during  the  second 

half  of  the  cyclotron  period  than  it  gains  during  the  first  half.  <>n  the 

whole,  its  energy  decreases  and  its  guiding  centei  shifts  to  the  left  (Fig. 

2b).  The  electron  energy  is  transferred  to  the  rf.  Electron  B  (Fig.  6c) 

is  initially  decelerated  by  the  if  field  and,  as  a  result,  its  tadius 

decreases.  Half  a  cyclotron  period  later,  electron  B  experiences  a  weakei 

acceleiating  electric  field.  Thus,  because  of  the  spatial  inhomogeneity  of 

the  rf  electric  field  cteated  by  the  lidges,  electron  B  loses  energy  to  the 

rf  during  a  cyclotron  period  and  its  guiding  center  drifts  to  the  right 

|Fig.  2 c  | .  Note  that  both  electrons  A  ai  d  B  yield  their  kinetic  energy  to 

the  rf.  Computet  simulation  shows  tha.  the  decelerated  electrons  can 

virtually  be  stopped,  leading  to  a  theoretical  efficiency  close  to  1<M) 
72  11 

percent  '  (  Experimental  results  may  he  moie  difficult  to  rntepret, 

however ‘  ) 

The  above  description  indicates  that  the  mechanism  of  energy  transfer 
in  peniotron  need  not  involve  beam  bunching  along  the  unpet tutbed  orbit  of 
an  electron  beam.  As  in  a  magnetron.  it  !■  the  transverse  migration  ul 
elections  into  regions  of  *  1 1  onge  i  t  I  t  t  e  1  if  created  fry  the  lidges  which 


leads  to  the  high  effic  iein 
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transverse  nodes.  A  nore  precise  definition  will  be  given  in  the  next 

section.  We  should  also  point  out  that  the  above  description  is  not  the 

only  possible  aechanism  responsible  foi  the  radiation  generated  in  a 
2 1 

peniotron  .  We  shall  say  more  about  this  at  the  end  of  this  section, 

where  we  treat  the  gyromagnet ron*0, The  gyromagnet ron  is  topologically 

identical  to  the  peniotron,  but  is  envisioned  to  operate  with  an  entirely 

different  mechanism. 
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Orbitron  is  similar  to  a  magnetron  in  that  an  electrostatic  field  is 

present,  and  that  the  rf  derives  most  of  its  power  from  the  electrostatic 

field.  A  model  of  orbitron  which  has  been  subject  to  some  theoretical 

study  is  shown  in  Fig.  7.  In  this  figure,  a  cylindrical  layer  of  electrons 

rotates  about  a  center  wire  which  is  positively  charged  with  respect  to  the 

outer  cylinder.  The  radial  electrostatic  field  provides  the  centripedal 

24 

acceleration  and  no  external  magnetic  field  is  present  .  In  some  sense, 

the  orbitron  may  he  considered  as  diametrically  opposite  to  the  peniotron, 

in  that  one  uses  only  an  electrostatic  field  for  electronic  motion  whereas 
the  other  uses  only  a  magnetic  field.  Moreover,  one  relies  on  complicated 
wall  structures  whereas  in  the  other,  a  smooth  wall  suffices,  in  which  case 
the  importance  of  tf  field  inhomogeneity  diminishes. 

Since  the  role  of  tf  field  inhomogeneity  is  minimal  in  an  orbitron, 
t  he*  only  way  for  the  beam  to  yield  its  energy  (kinetic  and  potential)  is 
via  a  hunching  mechanism  along  the  pat  tic  le  orbit  for  the  model  shown  in 
Fig.  1.  This  is  indeed  found  to  he  'tie  case.  Suppose  that  some  chaige 
perturbation  is  set  rip  on  the  beam.  A  test  election  moving  ahead  of  this 
charge  bunch  i-  a*  <  e  I  *■  i  a  '  c  d  .  if  out. .ltd  to  ,i  larger  radius.  This 

lead-  ’  ■  t  cl*-  i  *.  t  c  o  ;  t  In  i  *  •  t  ,i  t  j  on  ,  i  f  i  •  <pit  in  .  •-  l  in  e  the  i  egu  1  a  r  f  l  ecjuom  v 
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having  a  "negative  inertia"  -  it  falls  back  toward  the  bunch  upon 
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being  pushed  '  ’  .  Similarly,  a  test  electron  which  is  initially  behind 

the  charge  perturbation  is  decelerated,  it  moves  to  a  smaller  radius,  but 
its  angular  frequency  increases,  effectively  catching  up  with  the  bunch 
azimuthally.  Thus  there  is  strong  tendency  for  beam  bunching  in  the 
orbitron  configuration  shown  in  Fig.  7.  In  fact,  we  shall  show  in  the  next 

section  that  this  negative  mass  behavior  is  indeed  very  pronounced  - 

far  more  pronounced  than  that  exhibited  in  a  large  orbit  gyrotron  at  low 
beam  energy  ( < 1 0  KeV)  | See  Eq .  (27)  below). 

A  highly  bunched  beam  can  readily  yield  radiation.  In  so  doing,  the 
beam  contracts  radially  in  an  orbitron.  It  falls  to  a  region  of  lowei 
potential  energy.  It  can  easily  be  estimated  that  the  radiation  draws 
energy  primary  from  the  electrostatic  energy  like  in  a  magnetron. 

However,  as  the  beam  falls  toward  the  center  wiie,  its  angulai  frequency 
may  be  considerably  detuned  from  the  mode  of  opeiation,  thereby  limiting 
the  operating  efficiency.  Unless  some  form  of  phase  focusing  is 

introduced,  orbitrons  cannot  be  expected  to  yield  efficiencies  comparable 
to  magnetrons.  Thus,  the  model  shown  in  Fig.  7  has  the  peculiar  property 
that  it  yields  tremendous  small  signal  gain,  but  the  effi<  iency  may  be 
modes  t  . 

The  electron  beam  i  ri  a  large  oibit  gyrotron  |Fig.  Ha)  requires  an 

external  axial  magnetic  field  to  provide  the  rotation  There  is  no 

externally  i  mposed  IK  e 1 e>  t  i  i  <  t  i  e  1  d  ;  n  t  he  <  i  i  r  u  i  t  .  In  r  lie  s  i  mp  !  e  r 

geometry,  the  wall  corrugation  i-  absent  and  the  trail-;  let  of  energy  i 
through  charge  bunching  along  the  beam  orbit  jligs  Ha,  Bh|.  Similat  to 
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perturbation  in  the  electron  beam,  test  electron  A  in  front  of  the  bunch  is 


accelerated.  Its  rotational  frequency  |e |Bo/mQYo(A)  decreases  as  its 

relativistic  mass  factor  yq(A)  is  increased.  Azimuthally,  therefore,  A 

falls  back  to  the  bunch  as  a  result  of  this  relativistic  mass  effect. 

Similarly,  electron  B  is  decelerated,  yq(B)  decreases  but  its  angular 
frequency  <*>o(B)  =  |e|  B0/m0Y0(®)  increases.  Azimuthally,  B  speeds  up  and 
catches  up  with  the  bunch  |Fig.  8b].  The  negative  mass  instability  thus 
developed  was  shown  to  be  identical  to  the  cyclotron  maser  instability 
which  leads  to  radiation  in  gyrotions.  Howevet ,  it  will  be  shown  that,  at 
low  energy  beam  (<10KeV),  the  small  signal  gain  in  a  gyrofron  is 
considetably  smaller  than  the  otbirton.  On  the  othei  hand,  the  efficiency 
of  gyrotion  may  be  noticeably  highei  situe  the  magnetic  field  can  be 
detuned  so  that  w  w  is  sufficiently  small  to  yield  appreciable  gain,  but 
laige  enough  to  yield  substantial  efficiency  . 

Gy rot  tor  is  tuned  by  the  axial  magnetic  field.  The  magnetic  field 

requirement  can  be  substantially  reduced  if  the  tube  can  operate 

efficiently  at  a  high  eye  lotion  harmonic  .  A  promising  way  to  encourage 
harmonic  opera (ion  is  to  introduce  corrugation  at  the  waveguide  wall  |Fig. 
Hi  |  .  Tit  is  would  allow  effie  ient  operation  even  at  low  beam  energy.  Such 
an  1 n t  et  es  t i ng  possibility  wa  ■-  pi ed  r  c  t  ed  t  heoi e  t i ca 1 1 y  and  vei  i  f  i  ed 

expei  i  men  tall  v  when  ele<  t  i  on  beam-;  with  energy  at  tens  of  KeV  were 
employed'  to  ope i a  t e  at  high  ■  v  <  lotion  ha  i  mon i >  s .  The  resemblance  of  the 

magi  ie  1  i  on  on  t  e  i  .all,  t  oge  i  be  :  .  i  1 1 1  '  h.  g  v  i  o  »  i  on  ’in"  ban  1  sm  ,  I  €>ad  s  to  the 
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gyro-harmonictron,  or  harmonic  gyrotron,  have  later  been  used  for 

gyromagnetron.  But  none  of  the  latter  names  would  suggest  the  gyrotron 

mechanism  enhanced  by  a  magnetron  waveguide.  [Ref.  10  contains  a  lengthy 

description  of  the  historical  development  up  to  1982.) 

Comparing  Figs.  8c  and  6a,  one  immediately  sees  that  the  topologies  of 

the  gyro-magnetron  and  of  the  peniotron  are  identical!  Both  employ  annular 

beams,  ridged  waveguides,  and  only  an  external  magnetic  field  is  used  to 

provide  the  electron  rotation.  Even  the  beam  energy  and  frequency  range  in 

existing  experiments  are  similar.  The  immediate  question  is:  How 

important  is  the  transverse  peniotron  mode  in  the  operation  of 

gyromagnetron,  and  vice  versa.  Phrased  more  bluntly,  the  question  becomes: 

Is  the  peniotron  really  a  cyclotron  maser  in  disguise?  Or  does  a 

gyromagnetron  really  operate  according  to  what  is  envisioned  to  happen  in 

a  peniotron?  These  are  interesting  questions  which  have  generated 

29 

considerable  controversy,  and  a  resolution  is  not  immediate  .  Thus,  given 

the  magnetic  field,  beam  and  circuit  parameters,  one  cannot  even  predict 

21  29 

with  complete  confidence  in  which  mode  the  tube  should  operate  ’  ! 

If  the  simple  peniot ron/gyro  magnetron  configuration  can  cause 
uncertainties  in  our  understanding  of  their  operation  at  this  late  stage, 
just  due  to  the  corrugation  of  the  waveguide  circuit,  it  would  hardly  be 
surprising  that  the  situation  in  magnetron  would  be  much  worse  since  there 
is  not  even  a  consensus  in  tire  choice  of  a  suitable  equilibrium  model,  not 
to  mention  a  host  of  other  complex  features*’  and  the  highly  nonlinear 
behavior.  it  then  seems  desirable  to  restrict  here  to  a  tractable  thorny 
which  r  an  covet  'he  effects  of  <otiugated  structure,  curvature  and 
lelativtstic  mass  effects,  fur  an  eie  t  i  <>n  flow  undet  a  geneial  combination 
of  l>r  elec  t  i  i  c  field  and  magnetic  field,  with  the  efforts  of  t  lie  AC  and  DC 
self  field  included.  Tin ■  h  a  theory  will  be  presented  in  the  following 


sections  by  the  use  of  a  simple  model  (  Fig. 9  ).  There,  we  shall  isolate 
the  circuit  effects  and  the  intrinsic  beam  stability,  and  examine  their 
coupling. 


IV.  Stability  of  Laminar  Flows 


We  shall  now  give  a  simple  exposition  of  the  various  types  of 

instabilities  encountered  in  cross  field  devices.  These  instabilities  have 

been  invoked  to  explain  the  intense  noise  observed  and  the  initial  stages 

2 

of  spoke  formation  .  For  simplicity  of  mathematics  and  for  ready 
comparison  with  other  radiation  sources,  only  the  electron  sheet  model 
shown  in  Fig.  9  will  be  covered  in  some  detail.  Though  highly  simplified, 
this  sheet  model  retains  the  essential  features  of  the  physical  mechanism 
described  in  the  previous  sections.  It  also  demonstrates  the  profound 
effects  of  the  equilibrium  type  on  the  stability  of  an  electron  beam. 
These  effects  were  noticed  only  in  recent  years.  In  fact,  to  what  extent 
would  they  enter  in  the  case  of  a  thick  beam  remains  unanswered.  It  is 
highly  probable  that  the  effects  are  large,  since  different  waves  may  be 
excited  across  a  thick  shear  layer  and  they  may  interact.  We  shall  address 
some  of  these  issues  in  the  next  section. 

We  shall  begin  with  the  "longitudinal  mode".  The  real  part  of  the 
frequency  «  of  the  longitudinal  modes  satisfies 

«  =  /  <*>Q  (13a) 

where  ^  is  the  azimuthal  mode  number  and  coq  is  the  angular  frequency  of 
rotation  of  the  electron.  When  a  longitudinal  mode  is  present,  an  electron 
experiences  the  rf  field  at  an  almost  constant  phase  since  its  Doppler 
frequency  w-  /  =0,  and  there  would  be  a  strong  interaction  between  the 

electron  beam  and  the  rf.  This  is  particularly  true  if  w  is  also  close  to 
the  natural  mode  of  the  circuit.  The  longitudinal  modes  are  important  in 
gyrotron,  orbitron,  and  in  cross  field  devices. 

The  frequencies  of  the  "transverse  modes"  satisfy 

w  ■  /  Wq  -  ±  k  (13b) 

where  K  is  a  natural  frequency  of  oscillation  about  the  equilibrium  when  an 
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electron  orbit  is  perturbed,  [e.g.,  K  =  eBQ  /me  for  gyrotron,  peniotron, 

and  the  planar  magnetron,  but  K  =  <12  «  for  the  orbitron  model  (Fig.  7) 

where  an  electron  circulates  under  a  central  force  field,  generated  by  a 

positively  charged  thin  wire,  with  1/r  dependence].  The  peniotron  is 
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thought  to  operate  with  the  transverse  mode  *  ’  .  The  transverse  mode 

19 

is  also  an  important  one  in  cross  field  configurations 

Both  transverse  modes  and  longitudinal  modes  have  been  formulated 
27  28 

exactly  ’  for  a  cylindrical  laminar  flow  of  arbitrary  density  profile, 
including  DC  and  AC  space  charge  effects  for  a  general  combination  of 
radial  electric  field  EQ(r)  and  axial  field  BQ(r).  When  specialized  to  a 
thin  sheet  (Fig.  9],  a  simple  dispersion  relation  has  been  obtained.  The 
theory  presented  below  will  highlight  the  remarkable  features  which  are 
confirmed  by  a  direct  numerical  solution  of  the  exact  eigenvalue  problem. 

(A)  Longitudinal  Instability 

The  free  energy  which  drives  the  longitudinal  modes  unstable  may 

either  be  the  potential  energy  or  the  kinetic  energy.  The  potential  energy 

may  arise  from  the  DC  space  charge  of  the  electron  beam,  in  which  case  the 

30-36 

"diocotron  instability"  is  excited  ,  or  from  the  externally  imposed 

potential,  as  in  the  orbitron  model.  In  the  latter  case,  the  negative  mass 

instability  is  excited.  On  the  other  hand,  it  is  the  kinetic  energy  of  the 

electrons  which  is  converted  to  the  rf  in  gyrotrons.  As  we  have  explained, 

it  is  also  the  negative  mass  instability  which  is  responsible  for  the 

36 

radiation  generated  in  gyrotrons  .  The  transition  of  these  various 


instabilities,  and  their  relative  importance  in  different  regimes,  will  be 
covered  in  the  following  subsections. 


30-33  19 

The  diocotron  instability  ’  was  discovered  during  the  initial 
phase  of  magnetron  research.  It  was  also  one  of  the  first  plasma 
instabilities  ever  discovered.  The  free  energy  which  drives  this 
instability  resides  solely  with  the  electrostatic  energy  of  the  self  field 
of  the  electron  sheet.  There  are  several  (equivalent)  ways  to  view  this 
instability.  Two  of  them  are  given  below  for  the  model  shown  in  Fig.  10, 
which  is  a  further  simplification  from  the  model  in  Fig.  9. 

The  original  papers  attributed  the  diocotron  instability  to  a 

30 

hydrodynamic  instability  of  shear  flow  .  The  velocity  shear  is  due  to  the 

self  field  of  the  electron  sheet,  which  is  assumed  to  have  a  small,  but 

finite  thickness.  In  fact,  it  was  found  that  the  diocotron  instability  has 

the  same  growth  rate  as  the  Kelvin  Helmholtz  instability  of  a  vortex  sheet 

37 

in  classical  hydrodynamics  .  The  evolution  of  the  instabilities,  such  as 
the  curling  up  of  the  sheet,  is  also  similar. 

Associated  with  the  DC  self  field  of  the  electron  sheet  is  a  velocity 
shear  IFig.  10a]  in  the  equilibrium  ExB  drift.  The  velocity  differential 
due  to  self  field  is 


a 

o 


E  B  f 


o  o 


where  Es  is  the  total  change  of  the  self  electric  field  across  the  electron 
layer  of  thickness  x  and  constant  electron  density  pq.  Note  that  Av 
depends  only  on  the  surface  charge  density  oq,  but  is  independent  of  the 
thickness  x.  It  provides  a  good  approximation  for  a  general  density 

profile  if  X  is  sufficiently  small.  The  growth  rate  of  the  Kelvin 

37 

Helmholtz  instability  in  a  thin  vortex  sheet  is 
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where  k  is  the  wave  number  in  the  direction  of  flow.  Expression  (14)  is 
the  diocotron  growth  rate  of  a  thin  electron  sheet  in  the  long  wavelength 
limit  kx  <<  1.  Note  that  the  growth  rate  is  proportional  to  the  beam 
current . 

An  equally  transparent  way  to  see  the  physical  origin  of  the  diocotron 
growth  is  through  the  process  depicted  in  Figs.  10b, c.  Suppose  that  a 
sinusoidal  ripple  is  introduced  on  the  electron  sheet  [Fig.  10b). 
Electrons  in  positions  A,  C  experience  an  electrostatic  force  in  the  upward 
direction,  due  to  the  electrons  in  their  immediate  neighborhood.  Electrons 
in  positions  B,  D  experience  a  similar  force,  but  in  the  downward 
direction.  These  forces  produce  an  FxB  drift  to  the  left  for  electrons  at 
A,  C  but  to  the  right  for  electrons  at  B,  D.  This  leads  to  an  accumulation 
of  charges  between  B  C  and  a  deficiency  of  charges  between  A  B  [Fig.  10c). 
The  electric  field  associated  with  these  charge  perturbations,  as  shown  in 
Fig.  10c,  produces  an  ExB  drift  of  electrons  which  reinforces  the  original 
ripple,  and  the  instability  grows  as  a  result. 

The  diocotron  effect  was  considered  to  be  a  major  cause  for  the 

intense  noise  and  fluctuation  observed  in  cross  field  devices  and  high 

power  diodes.  It  was  also  thought  to  contribute  to  a  significant  increase 

2  3 

of  the  gain  in  the  theories  of  cross  field  amplifiers  ’  .  There  are  many 
refinements  of  the  diocotron  instabilities,  such  as  the  effects  of  finite 


shall  stick  to  the  thin  layer  model  [Fig.  9]  so  as  to  bring  out  the  crucial 


dependence  of  stability  on  the  equilibrium  type.  For  example,  we  shall 

show  below  that  the  diocotron  instability  may  be  stabilized  if  the 

curvature  effects  are  taken  into  account, 
i 

I 

! 

(A2)  General  longitudinal  instability  without  DC  space  charge  effects 
[negative  mass,  cyclotron  maser,  orbitron  modes,  etc.]. 

The  curvature  and  relativistic  effects  have  been  ignored  in  the  above 
description  of  the  diocotron  instabilities.  In  practice,  in  many  radiation 
sources  the  curvature  effects  of  the  electron  orbits  are  crucial.  The 
relativistic  effect  may  already  be  important  for  an  electron  beam  with 

i 

energy  as  low  as  5  KeV,  as  shown  in  the  original  experiments  on  cyclotron 
38 

maser  .  To  examine  all  of  these,  we  now  pretend  in  this  subsection  that 
the  electron  layer  is  neutralized  to  rule  out  the  diocotron  instability. 
The  questions  become:  Does  it  make  a  difference  in  the  stability  if  the 
equilibrium  rotation  is  supported  by  a  radial  electric  field  alone  (like  an 
orbitron),  or  by  an  axial  magnetic  field  alone  (like  a  large  orbit 
gyrotron)  or  by  a  combination  of  both  (like  a  smooth  bore  magnetron)?  How 
and  in  what  way  would  the  relativistic  effects  enter?  All  of  these 
questions  can  now  be  answered  with  the  following  simple  study. 

Consider  a  thin,  neutralized  electron  layer  rotating  concentrically  at 
velocity  vq  =  9  v  (r)  =  9rwo(r)  inside  some  waveguide  structure 
[Fig.  9).  The  E-layer  has  a  small  thickness  t,  and  uniform  charge  density 
extending  from  r=r^  to  r=r£  mean  radius  R.  For  the  time  being,  we 

shall  leave  unspecified  the  relative  strength  of  the  radial  electric  field 
Eq  and  the  vertical  magnetic  field  which  are  needed  to  provide  the 

circular  motion  of  the  beam.  Thus,  vq  is  governed  by 
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It  is  convenient  to  introduce  a  dimensionless  quantity  h  to  characterize 
the  strength  of  the  electric  field  of  the  equilibrium: 


h  = 


-erE 


3  2* 
m  y  v 
oo  o 


(16) 


2 

Note  that  YQh  is  equal  to  the  ratio  of  the  electric  force  to  the 
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centrifugal  force  in  equilibrium  .  Using  this  definition  and  the 
equilibrium  condition  (15),  one  may  characterize  the  equilibrium  type 
according  to  the  value  of  h  [Fig.  11]: 


(i)  h=0  corresponds  to  the  large  orbit  gyrotron  [Fig.  8],  where  the 
equilibrium  rotation  is  supported  by  an  axial  magnetic  field  alone  (Eq=0). 

2 

(ii)  h=l/yo  corresponds  to  the  orbitron  model  [Fig.  7],  in  which  the 
rotation  is  supported  solely  by  a  radial  electric  field.  [Bq=0,  cf.  Eqs. 
(15),  (16)]. 

2 

(iii)  h»l/ro  corresponds  to  an  inverted  magnetron,  with  the  cathode  at  the 
outer  conductor  and  the  anode  at  the  inner  conductor.  The  rotation  is 
approximately  given  by  the  ExB  drift  (Centrifugal  force  is  small). 

2 

(iv)  h<<-l/yo  corresponds  to  a  conventional  magnetron,  with  the  cathode  at 
the  inner  radius  and  the  anode  at  the  outer  radius.  Again,  the  rotation  is 
approximately  given  by  the  ExB  drift  and  the  centrifugal  force  is  small 
compared  with  either  the  electric  or  Lorentz  force  in  equilibrium. 


(v)  The  planar  limit  is  recovered  formally  as  r  ■*  «  (fixing  Eq,  vq).  That 
is,  |h  |  -♦  “corresponds  to  the  planar  limit  [cf.  Eqs.  (15),  (16)1. 

These  various  special  values  are  labelled  on  the  h  axis  in  Fig.  11  and 
a  simple  dispersion  relation  is  derived  below  for  general  values  of  h. 
Since  we  are  examining  only  the  longitudinal  modes,  the  beam  interacts  with 
an  rf  field  mainly  through  the  azimuthal  component  of  the  electric  field 
(Ejq)  which  the  thin  beam  experiences.  Conservation  of  energy  gives 

ev  E,  Q  =  de/dt  (17) 

o  i  y 

where  t  is  the  total  energy  (kinetic  and  potential)  of  the  electron  beam. 
Upon  using  the  chain  rule,  we  express 

de/dt  =  (d0/dt)/(d9/de)  =  (h/R)/(d0/de) 
in  terms  of  the  azimuthal  displacement  h  of  an  electron  from  its 
unperturbed  position.  Thus,  (17)  becomes,  upon  linearization, 


D  =  eRv 


rlw  . 

oEi0tar)  5 


eE10/Meff 


where  the  equilibrium  value  0=wq  is  expressed  as  a  function  of  the  particle 

energy.  In  analogy  with  the  force  lav  "F=ma",  we  define  in  (18)  an 

effective  mass  M  cc  =  (Rv  dw  /de)  It  is  not  difficult  to  show  from  (15) 
eff  o  o 

and  (16)  that 

»e££  ‘ 

o 


where  0  =  v  /c. 

o  o 


Equations  (18)  and  (19)  are  very  interesting.  They  govern  the 
longitudinal  (azimuthal)  dynamics  of  the  electron  beam.  The  crucial 
properties  of  gyrotron,  orbitron,  gyromagnetron ,  smooth  bore  magnetron  and 


its  inverted  configuration,  or  even  cyclic  particle  accelerators  such  as 
betatrons  are  contained  in  this  effective  mass.  Relativistic  effects, 
curvature  effects,  and  the  magnitude  and  sign  of  h  determine  this  effective 
inertia,  which  in  turn  determines  the  response  of  an  electron  to  an 

accelerating  or  decelerating  rf  electric  field.  For  example,  is 

2  2  2 

negative  for  h>-6  / 2,  is  positive  when  h<-g  / 2 ,  is  infinite  at  h=-0  /2,  and 
o  o  o 

2 

has  a  local  maximum  when  h=l/vo  (Fig-  1 2 | .  All  of  these  are  interesting 
properties  and  will  be  considered  shortly.  Note  for  the  time  being  that 
the  main  dynamical  properties  associated  with  the  equilibrium  type  for  the 
longitudinal  modes  are  already  contained  in  (18),  (19). 

For  a  perturbation  proportional  exp|  j  wt  -  jj[  9j ,  d/dt  stands  for 
j(«-(wo)  where  is  the  azimuthal  mode  number.  Associated  with  the 

azimuthal  displacement  n  is  a  surface  charge  density  perturbation  ,  given 
by 
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(20) 


in  which  is  understood  to  be  the  value  evaluated  at  r=R.  In  writing 
the  last  expression,  we  have  used  (18).  Equation  (20)  expresses  the  charge 
density  perturbation  in  the  beam  in  response  to  some  imposed  azimuthal 
electric  field.  Thus  far,  the  surrounding  electromagnetic  structure  has 
not  entered  into  consideration. 

To  complete  the  analysis,  we  calculate  what  kind  of  electromagnetic 
waves  would  be  excited  if  there  is  an  RF  surface  charge  located  at 
r-R.  The  piesence  of  a  surface  cliatge  ,  by  Gauss  law,  produces  a 


discontinuity  in  the  tadial  component  of  the  RF  electric  field  E.  .  Thtis 


where,  and  in  what  follows,  the  supei script  .(  )  is  used  to  designate  the 

position  just  outside  (inside)  of  the  tadius  undei  consideration.  We  next 


express  *  in  teres  of 
at  the  outer  and  inner 
defined  by 


edges 


by  introducing  the  wave  admittances  b  and  b 
of  the  beam.  Specifically,  b  and  b  ate 
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Since  the  tangential  electric  fie  Id  (E10  )  is 
across  the  electron  layer,  we  have  E^  *  = 

give 


i 

i 

i 
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a 

e 


1 
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•  b  ) E j  q. 


approximately  continuous 
=  E^q  and  hence  (21)  (21) 


(  24  ) 


Whereas  Eq .  (20)  expresses  the  dynamical  l espouse  of  the  beam  to  an 
rf  electric  field  without  reference  to  the  circuit,  Eq .  (24)  gives  the 
circuit  response  to  some  beam  excitation  (Oj),  irrespective  of  how 
is  produced  dynamically.  For  .elf  consistency,  we  obtain  from  (2<>),  ( 2* ) 
the  dispersion  relationship 
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In  (26),  N-2nRoo/e  is  the  charge  number  density  per  unit  axial  length,  \»  is 


the  dimensionless  'Budkei  parameter'  measuring  the  beam  density  (cut  tent) 


I  VsNe^/^llt  m  e‘  in  MK  S  units,  it  equals  to  e^N-m  '  m  (<;S  un  1  f  s .  I 
1  o  O  n  1 


The  dispeision  relationship  (.’b)  governs  t  lie  longitudinal  stability 


of  a  model  election  sheet  in  otbitions,  gyiotions,  gy i omagnet t on ,  and 

7, 


(smooth  bote  )  magnetrons.  The  effects  of  the  waveguide  wall  enter  only  in 
the  normalized  wave  admittance  h  ♦  b  in  t.’S)  and  the  dynamics 

(  ecru  i  1  l  hi  i  urn  t  vpe  )  enters  m  the  factor  (  .  ,’h  )  (  1  •  >  ‘  )  ,  wli  l  c  h  is  also  t  fie 

o  o 

ma  s  s  factor  in  F  t  j  .  (  1  '* )  .  The  wave  a  dm  i  t  •  am  *  l>  and  h  is  a  function  of 

geometty  and  of  frequency  ui,  and  t<>!  t  he  longitudinal  modes  we  ate 

studying,  ui=  /  vo  The  mathematical  expressions  of  h  rod  b  foi  various 

*  ^ 

cases  aie  given  in  Appendix  B  of  thernin  and  Ian'  The  majoi  points  ate 

as  follows.  Uhen  t  lie  waveguide  wall'  at*-  lossles?.  h  •  h  aie  teal.  It  can 

eithet  be  positive  u  negative.  'When  .'  i  •  po>i»i.e  (negative),  the  beam 

1 1 

l  said  t  o  expe  t  t  eii'  e  a  ■  a  pa  l  t  i  .  e  (  i  ml  m  •  i  .  c  )  ■-  t  ;  in  t  n  i  e  Fo  i  f  i  equenc  y 


ui  '  lose  to  cut  off  f  i  equenc  ies  (  )  of  the  vacuum  waveguide  mode. 


)  •  h  -» <  t*>  cw  ),  and  the  growth  iate  becomes  laige,  tegardless  of  the  s  i  jjn  of 


H  ...  Under  f  b  i  •;  condi  iron  (synchronous  condition),  (?'.)>  becomes  a  cubic 
e  f  f 


polynomial  of  (  w  Q  ^  ).  leading  to  a  i  a  t  of  giowtli  (gain)  proportional 


to  i  ,  which  i  <.  '  ha  i  a>  t  t>i  i  t  i  <  •  ■  t  gain  in  1  ong  i  ’  'id  i  no  1  i  n  t  ei  a<  t  l  ons  ,  as 

-01 ,  ih.c  1  ,  H 


n  g,  i  o  t  i  on  '•  '  ’  Hit"  c-  t.i'c  nc  tit  on  '  t  ic  v  a  i  a  dm  l  t  t  an-  e  s  ma  v  he 

«>v  f  ended  to  *  he  a  s  e  1  -fie  i  c  » fie  -a  .  ego :  •(»•  -  a  i  i  he.  c*m»  '  m  l  ugated,  as  in  a 
g',1 1  c  imagne  t  r  c >ii 

* ■<  i  ■  'h  r i . » ”•  i  .  we  'hall  a  ■:  ume 


F  o  t  '  ><n '  <  n  i  eri' t  -if  >f  . 
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(a)  Vhcn  h»0,  Eq .  (2b)  indicates  that  an  instability  exists  as  a  tesult 


of  the  relativistic  effect.  This  is  the  veil  known  negative  aass 
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instability  *  ,  whose  gtowth  tale  is  given  by  w.  (£t/R)0o  /(b  »b  ) 

which  vanishes  in  the  nonrelat i vist ic  limit  6  *o.  This  instability  was 

o 

predicted;  and  it  places  a  limit  on  the  beam  rurient  in  cyclic  accelerator. 
It  turns  out  to  be  identical  to  the  cy<  lotion  maset  instability4^,  which  is 
tesponsible  foi  the  ladiation  geneiated  in  gytotions. 


'j 

(b)  The  negative  mass  gtowth  rate  would  lie  present  as  long  as  h>  3^  >2. 

7 

It  is  maximized  with  tespect  to  h  when  hl/>",  as  teadilv  demonstrated  f  i  om 

o 

Eq.  (25).  This  case  cotiesponds  to  the  oibitron  configuration  |Figs. 

1 1 . 1 2 | .  In  other  words,  fot  a  given  totational  enetgy  and  a  given 

geometry,  the  negative  mass  effect  is  most  pronounced  when  the  equilibrium 

2  7,28 

is  supported  by  a  tadial  elec* tic  field  alone  as  in  an  orbit  ion  model 
By  setting  h  =  0  and  h=l/y^  in  (2b),  we  obtain  the  comparison  of  the  small 
signal  growth  between  a  large  orbit  gyrotton  and  an  or  bit  ion 


w.  (orbit r on)  ^ 

<*»•  (gvrot  ton)  6 

i  o 


(21) 


Since  0  equals  to  0.044,  o . 1 4 .  0.17,  o . bb ,  0.78  when  the  rotational 

o 

enetgy  is  IKeV,  SKeV,  40KeV,  lOOKeV,  and  U  lOKeV ,  respectively,  one 

sees  that  the  small  signal  gain  for  an  electron  layet  in  an  oibitron 

con  figuration  'an  in  principle  he  r  a  t  het  high  in  comparison  with  the  other 

radiation  sources,  especially  at  low  beam  energy-  A  potential  use  of  this 

4  1 

enhanced  charge  bunching  in  klystrons  has  been  suggested  recently 


(c)  The  dispet  sion  i  elation  (  '  > )  suggests  that  the  negative  mass 
instability  is  suppressed  if 

h  <  e‘  <.’«) 

<1 

That  is,  the  negative  mass  instability  may  he  stabilized  by  a  negative 
tadial  (X  eletiiit  field  of  a  suitable  magnitude.  In  teims  of  an 
exteinal  potential  V  imposed  between  the  timet  < ondin  to:  at  i  a  and  the 
oil  t  ei  <  ondin  t  <u  a  t  i  b  |  See  big  1  *  tie  1  ov  |  ,  ’  he  stability  <  ond  l  t  i  on  (  >H  ) 
i  ead s 
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Physically,  when  h=  0^/2,  the  effective  mass  of  a  totaling  electron  is 
infinite  |cf.  Pig-  12|.  The  beam  is  very  rigid  azimuthally  and  is 
incapable  of  transferring  its  rotational  energy  to  the  resistive  wail, 
which  is  the  physical  mechanism  fm  the  excitation  of  the  resistive 
i ns  t  abi  1  i  t  . 

(e)  The  negative  mass  instability  should  disappear  in  the  planar  geometry 
limit,  as  expected  intuitively.  This  is  also  reflected  in  the  dispersion 
relation  (25).  In  this  planar  limit  R  ♦CT\  /  ♦  »>  ,  but  jt  /  R,  T,  E  remain 

finite.  Then  h  ♦*  by  (lb)  and  the  r  ight  hand  member  of  (25)  tends  to  zero. 
Vhat  remains  is  then  the  diocottort  instability  which  is  not  included  in 
(25)  since  the  sel.  field  has  been  explicitly  ignored  in  the  derivation 
of  (25). 

Viewed  slightly  differently,  we  may  regard  (25)  as  the  lowest  order 

dispersion  relation  when  we  expand  the  growth  rate  in  terms  of  the  small 

parameter  (r/R).  The  diocotion  growth  rate  (14)  is  the  residual 

lb 

instability  when  the  rmvatute  effect  is  absent  This,  in  fact,  is  found 

to  he  the  rase  in  an  exact  formulation  of  the  eigenmodes.  This  will  be 
discussed  next,  and  examples  will  be  given  to  demonstrate  the  validity  of 
'be  growtfi  rate  formulas  (25)  thi-nigh  a  direct  numerical  integration  of  the 
exa-  '  eigenvalue  problem. 

<  A  '•  )  An  ex  a  i  formulation  nod  >  h»  «  i  1  *  i  .  a '  l  on  of  r  lie  simplified 
d  l  pe  i  ■  i  on  r  e  1  a  t  i  oil 
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combination  of  radial  electric  field  and  axial  magnetic  field  with  self 


consistent  density  and  velocity  profiles.  It  is  fully  relativistic  and 
fully  electromagnetic;  it  includes  both  DC  and  AC  space  charge  effects. 
For  TE  modes  with  dependence  exp(  jtot-j  £  0) ,  the  eigenvalue  equation  for 
♦*rEjg  is  of  the  form 

(Jr  (A  a?)  * B*  -  °-  <30> 

where  A  and  B  are  rather  complicated  functions  of  to,  £,  and  r.  This 

eigenvalue  equation  governs  both  the  longitudinal  modes  and  the  transverse 

modes.  These  modes  manifest  themselves  as  singularities  of  the  equation. 

The  complex  eigenf requencies  u*=to(/)  for  both  types  of  modes  are  obtained  by 

matching  the  solutions  of  $  at  the  beam  edges  to  those  of  the  vacuum 

solution  exterior  to  the  beam  [or  $=0  at  the  conducting  boundary  with  which 

the  electron  sheath  happens  to  be  in  contact.) 

The  eigenvalue  equation  (30)  has  been  solved  both  analytically  and 

47 

numerically  for  a  thin  beam.  In  the  analytic  procedure  ,  we  expand  all 

quantities  in  powers  of  (t/R).  To  two  orders  in  (x/R),  the  growth  rate 

.  27,28 
is  given  by 
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All  quantities  in  this  expression  has  been  defined  in  the  previous  section. 

Note  that  the  first  term  of  (31)  is  the  same  as  the  rhs  of  (25)  except  for 

39  0 

the  generalization  of  h  in  (31).  In  the  planar  limit  R-»®,  /  -**>, 

i/R-»k  =  finite,  then  h-**>,  and  the  first  term  of  (31)  disappears.  Using  the 


last  term  of  (31),  and  noting  that  only  the  term  involving  remains  in 
this  planar  limit,  ve  obtain 

2  ..2  2.,  4,  4.  2 
1  '  o'  p  c 

This  equation  is  just  the  relativistic  version  of  the  diocotron  instability 
for  a  sheet  beam  [cf.  Eq.  (14),  also  Ref.  (33)  J .  Thus,  the  diocotron 
instability  is  in  fact  recovered,  and  the  negative  mass  instability 
removed,  in  the  planar  limit.  In  other  words,  Eq .  (31)  includes  both  the 
DC  and  AC  space  charge  effect. 

For  a  thin  beam,  the  second  term  in  (31)  is  using  negligible  compared 
with  the  preceding  term.  Thus,  stabilization  of  the  negative  mass 
instability  by  the  use  of  the  criterion  (28)  would  also  imply  stabilization 
of  the  diocotron  instability  for  a  thin  beam.  Figure  13  shows  the  growth 
rates  according  to  the  dispersion  relation  (31)  in  two  examples.  The 
interesting  features  are  noted  there:  (a)  maximum  small  signal  growth  for 


orbitron  configuration  ( h= 1 / yq ) ;  (b)  stabilization  of  negative  mass  and 

2 

diocotron  instability  when  h<-8Q/2;  (c)  persistence  of  small  signal  growth 

at  low  beam  energy  for  gyrotron  (h=0).  This  growth  is  due  more  to  the 
diocotron  effect  than  to  the  negative  mass  effect.  [See  Ref.  36  for  more 
detail . ) 


The  validity  of  the  analytic  theory  has  been  confirmed  by  comparing 
the  growth  rates  with  those  obtained  from  direct  numerical  integration  of 
the  governing  differential  Eq .  (30).  Shown  in  Figs.  14,15  are  some 
examples  of  this  comparison.  The  solid  curves  represent  the  eigenvalues 
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according  to  numerical  integration  of  (30)  and  the  dashed  curves  to  (31). 

Note  that  the  stability  condition  and  the  peak  growth  rate  predicted  for 
orbitrons,  according  to  (25),  are  confirmed  in  Fig.  14.  The  agreement  in 
the  dependence  on  beam  current  (Fig.  15)  is  remarkable.  Other  checks  on 


geometry,  on  mode  number,  on  yq>  on  x/R,  etc.,  have  been  performed 
They  all  confirm  the  validity  of  the  analytic  dispersion  relation  (31). 
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B  Transverse  Instability 

The  longitudinal  modes  described  in  the  previous  section  rely  heavily 
on  the  spatial  inhomogeneity  of  the  unperturbed  motion  {e.g.,  velocity 
shear  due  to  space  charge  in  the  diocotron  instability,  and  differential 
rotation  in  the  negative  mass  instability  as  reflected  in  the  effective 
mass  factor  d(oQ/de=(d«o/dr )/(de/dr ) j .  The  spatial  inhomogeneity  of  the  rf 

field  plays  a  secondary  role.  Near  synchronism,  the  growth  rates  of  these 

1/3 

modes  are  proportional  to  I  .  In  contrast,  the  gain  mechanism  in  both 

cross  field  configuration  and  peniotron,  the  spatial  inhomogeneity  of  the 

rf  field  plays  a  crucial  role.  This  situation  arises  since  the  transverse 

drift  of  the  "favorable"  electrons,  acted  on  by  the  rf  fields, 

progressively  populate  a  region  with  stronger  rf  fields,  whereas  the 

"unfavorable"  electrons  drift  to  a  region  of  weaker  rf  field.  The  small 

1/2 

signal  gain  for  the  transverse  mode  is  proportional  to  I 

For  the  present  sheet  beam  model  [Fig.  9],  the  transverse  modes  are 
characterized  by 


/  5  \  ,,  2,2 . 1/2 
WWo>  =  ±  (1  +  Y0h  )  <oo. 


(33) 


2  2  1/2 

One  can  readily  show  that  the  frequency  (1  +  yoIi  )  coq  represents  the  radial 
oscillation  of  an  electron  about  its  guiding  center  in  a  frame  co-moving 
with  the  wave,  tinder  the  combined  action  of  E  and  B  .  For  example,  in  the 
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nonrelativistic  limit  of  orbitron  configuration  (h=l/YQ-»l),  this  frequency 
becomes  /2  wq,  which  is  the  frequency  of  radial  oscillation,  about  the 
guiding  center,  of  a  particle  under  a  central  force  field  whose  strength 

behaves  as  1/r  [Fig.  7].  In  another  limit  h=0,  (33)  becomes  w- 

n  ...  „  2,3,11,19,28-33,40,41 

j([«c=±toc  which  represents  gyro-resonance 

As  mentioned  earlier,  the  transverse  modes  are  also  governed  by  the 

same  differential  equation  (30).  Using  the  same  analytical  procedure  which 

leads  to  the  dispersion  relation  (31),  we  obtain  the  dispersion  relation 
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for  the  transverse  modes  for  a  thin  electron  layer  : 


<►/„., ^  (i2$  (l-  g  [£  -q)  £♦,)].  <34, 


Here 


q  a  -  <o 


,(1  +H)/[  (ft  -  (£«/  <0Q)  («  -  yt»0)  1  - 


(35) 


This  is  the  lowest  order  dispersion  relationship  in  the  expansion  parameter 
t/R.  It  is  valid  for  general  combination  of  Eq  and  Bq;  and,  since  in  this 
case,  DC  space  charge  is  not  critical,  it  gives  the  correct  result  in  the 
limit  of  planar  geometry. 

This  dispersion  relation  is  rather  complicated.  The  wave  admittance 
b+  and  b  reflects  the  corrugation  of  the  walls  and  the  rf  field 
inhomogeneities  such  corrugation  produces^' ^ '  32 , 41 .  without  going  into 
the  details,  we  may  make  the  following  remarks. 


(a)  In  the  absence  of  wall  loss,  small  signal  growth  can  occur  only  when 
the  transverse  mode  is  synchronous  with  the  circuit  mode,  for  which 


b  +b  =  0.  The  growth  rate  under  such  a  synchronous  condition  is 


.1/2 


proportional  to  I1'1'.  This  small  signal  growth  is  typically  weaker  than 
that  associated  with  the  longitudinal  modes. 
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(b)  Since  synchronism  is  always  required  for  small  signal  growth  the 
relative  importance  among  the  centrifugal  force,  Lorentz  force,  or 
electrostatic  force  in  supporting  the  equilibrium  seems  to  be  less 
important.  That  is,  the  transverse  mode  is  perhaps  less  sensitive  to  the 
equilibrium  type.  But  the  resistive  instabilities  are  equally  serious  for 
both  transverse  modes  and  longitudinal  modes. 

(c)  Conventional  theories  of  transverse  modes  also  express  the  gains  in 

terms  of  the  rf  electric  fields  (and  their  gradients)  at  the  beam 
3  A  41 

location^  ’  .  Equation  (33)  is  a  generalization,  including  a  possible 

radial  electric  field  in  a  curved  geometry. 


(d)  Numerical 

simulations 

indicated 

that  the  transverse 

modes , 

in 

principle,  can 

be  a  very 

efficient 

mode  to  operate 

at  as 

in 

22  29 

peniotrons  ’  .  The  small  signal  theory  only  gives  the  initial  stage  of 
the  development.  The  effects  of  the  spatial  inhomogeneities  of  the  rf 
field  are  incorporated  in  b+  and  b  ,  which  are  complicated  expressions 
themselves  when  ridges  are  present  in  the  waveguide  walls.  Much  work, 
remains  in  the  study  of  the  transverse  modes,  especially  their  possible 
coupling  with  the  longitudinal  modes  when  the  sheath  thickness  is  no  longer 
small . 
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V. 


Remarks 


A  rudimentary  theory  is  given  for  a  class  of  radiation  sources  in 
which  the  basic  electron  flow  may  be  regarded  as  laminar.  Emphases  have 
been  placed  on  the  dynamical  dependence  of  the  equilibrium  type  and  on  the 
role  of  circuit  structure.  In  general,  the  spatial  inhomogeneities  of  the 
rf  fields  are  important  when  the  transverse  migration  of  the  electrons  is 
essential,  as  in  cross-field  devices  and  in  peniotron.  On  the  other  hand, 
they  are  less  important  if  the  radiation  draws  upon  the  charge  bunching 
along  the  particle  orbits,  as  in  the  (large  orbit)  gyrotron  and  orbitron. 
In  the  latter  cases,  the  spatial  inhomogeneities  of  the  unperturbed  motion 
(i.e.,  velocity  shear)  are  far  more  important  than  those  of  the  rf  fields. 
Regardless  of  the  dynamical  properties,  all  of  these  radiation  sources 
always  operate  at  a  frequency  close  to  a  natural  frequency  of  the  waveguide 
circui t . 

Several  interesting  properties  were  predicted  from  a  study  of  the 
longitudinal  modes  on  a  thin  beam.  For  example,  the  orbitron  configuration 
offers  perhaps  the  highest  small  signal  gain,  especially  at  low  beam 
energy.  Away  from  synchronism,  both  the  negative  mass  and  diocotron 
instabilities  may  be  stabilized  by  a  negative  radial  electric  field  of 
sufficient  magnitude.  This  arises  as  the  electron's  effective  mass  is 
changed  from  negative  to  positive.  [We  should  stress  that  instability  will 
occur  if  the  circuit  mode  is  synchronized  with  the  beam  mode,  even  if  the 
beam  has  a  "positive  mass"  behavior.]  The  most  stable  configuration 
against  synchronous  excitation  and  against  even  the  resistive  wall 
instability  is  when  the  beam  is  very  "rigid".  This  occurs  when  the 

effective  mass  is  infinite  -  i.e.,  when  the  radial  electric  field  is 

adjusted  so  that 

h  -  -  g2/2  . 
o 
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All  of  above  predictions  are  yet  to  be  confirmed  in  controlled 
experiments,  such  as  the  one  proposed  in  Fig.  16.  This  is  essentially  a 
two  cavity  klystron  with  a  bent  drift  tube.  The  orbit  is  bent  either  by  a 


magnetic  field  Bq  or  by  an  electrostatic  field,  or  by  both.  The  polarity 

and  the  magnitude  of  the  externally  imposed  voltage  (V)  and  the  external 

magnetic  field  Bq  may  be  adjusted  to  correspond  to  various  values  of  h 

[Fig.  11].  The  response  of  the  electron  beam  may  be  monitored  at  the 

output  cavity,  after  an  external  rf  signal  is  impressed  upon  the  beam  at 

the  input  cavity.  Such  an  experiment  may  be  carried  out  with  an  electron 

beam  of  energy  <  10  KeV,  Bq  <  100G,  Eq  <  5  KeV/cm,  and  beam  current  <  0.1A. 

It  also  serves  as  a  proof-of-principle  experiment  for  the  "super-bunched" 
43 

klystron  ,  as  veil  as  a  controlled  experiment  on  the  orbitron  mechanism 
discussed  in  this  Chapter. 

In  the  case  of  a  thin  beam,  the  longitudinal  mode  is  decoupled  from 

the  transverse  mode,  and  each  mode  may  be  treated  separately.  This  is  not 

the  case  when  the  beam  is  thick.  Within  a  thick  beam,  because  of  the 

velocity  shear,  both  Eqs.  (13a)  and  (13b)  may  be  satisfied  for  electrons  at 

different  radii.  Thus,  the  longitudinal  mode  of  the  electrons  at  one  layer 

may  interact  with  the  transverse  mode  of  the  electrons  at  a  different 

layer.  In  fact,  this  was  noticed  in  the  original  stability  theory  of  the 

19 

Brillouin  flow,  and  Buneman  in  effect  already  interpreted  the  instability 

as  due  to  an  exchange  of  the  negative  energy  wave  between  those  electrons 

satisfying  (13a)  with  the  positive  energy  waves  of  electrons  at  a  different 

layer  at  which  (13b)  is  satisfied!  A  rich  amount  of  interesting  physics 

buried  in  Eq.  (30)  remains  to  be  explored. 

There  are  deficiencies  in  the  laminar  flow  model.  The  lack  of  Larmor 

motion  in  such  a  model  was  only  recently  corrected  in  the  treatment  of 

35 

diocotron  instability.  Kleva,  Ott  and  Manheimer  showed  that  the 
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diocotron  instability  cannot  be  completely  removed  from  an  electron  sheath, 

regardless  of  the  relative  size  between  the  sheath  thickness  and  the  Larmor 

radius.  [This  prevalence  of  the  diocotron  instability  has  recently  been 

conjectured  to  occur  within  the  Debye  sheaths  of  a  magnetized  plasma  and 
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may  give  rise  to  the  phenomenum  of  Bohm  diffusion  . J 

Dynamically,  perhaps  the  most  serious  deficiency  in  existing 
magnetron  theories  is  the  failure  to  take  into  account  the  periodic 
variation  in  the  DC  electric  field  resulting  from  the  wall  corrugation. 
Such  a  periodic  variation  should  also  be  introduced  in  the  unperturbed 
state  of  the  electron  flow.  If  one  recalls  the  importance  of  the  periodic 
motion  in  the  mean  flow  of  the  electrons  in  free  electron  lasers,  one 
cannot  fail  to  worry  about  the  possible  omission  of  significant  physics 
when  some  built-in  "wiggling  motions"  have  entirely  been  ignored  in 
specifying  the  unperturbed  state  of  a  magnetron  in  virtually  all  existing 


theories  . 
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Fig.  1  (a).  Laminar  flow  model  of  magnetron. 


(b).  Further  idealization. 


VOLTAGE  (V) 


Fig.  2  The  Hull  cut-off  voltage  and  the  Buneman-Har t ree  threshold  voltage 
The  latter  curve  is  tangent  to  the  former. 


^-CATHODE 


THEORETICAL  BRILLOUIN  LAYER 

Fig.  A  Same  as  in  Fig.  3,  but  at  a  later  stage  when  oscillation  is  fully 

developed.  The  magnetron  voltage  is  reduced,  leading  to  a  narrower 
Brillouin  sheath  as  shown  in  this  figure. 


/y 


(a) 


Fig.  6  (a)  A  peniotron  model. 

(b)(c)  Action  of  the  rf  electric  field  on  electron  A  which  is 

initially  in  an  accelerating  phase  and  on  electron  B  which 
is  initially  in  a  decelerating  phase.  The  dotted  circles 
denote  the  unperturbed  orbits. 
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7  An  orbitron  model.  The  electron  layer  rotates  about  the  center 
conductor  under  a  radial  electric  field.  This  model  is  used 


in  Refs.  (12,  27,  28). 
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Fig.  8  (a)  A  gyrotron  model. 

(b)  Negative  mass  effect. 

(c)  A  gyromagnetron  model.  A  magnetron  waveguide  is  used  to 
encourage  harmonic  generation. 
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NON-NEUTRALIZED  ELECTRON  SHEET 


.10  Physical  origin  of  the  diocotron  instabilities,  viewed  in 
the  rest  frame  of  the  mean  flow: 

(a)  velocity  shear  generated  by  the  electron  self  field. 

(b) (c)  ripple  on  the  sheet  and  its  reinforcement. 
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Fig. 11  Correspondence  between  the  values  of  the  normalized  electric  field 
h  and  the  various  types  of  equilibrium. 
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Fig. 13  Growth  rates  of  longitudinal  modes  in  a  thin  electron  layer  as  a 
function  of  the  equilibrium  type  ( h ) ,  at  both  non-relativistic 
energy  (top)  and  relativistic  energy  (bo t tom)  They  are  calculated 
from  (31)  with  [  =1,  t/R=0.016,  b+=5,  b_=3,  and  beam  kinetic 
energy  and  density  as  specified  in  the  figure. 


Fig. 14  Normalized  growth  rate  for  longitudinal  modes  as  a  function  of  h 


for  a  test  case  with  a  =  0.6m,  b=2.6m,  R  =  1.0m,  t/R  =  0.02, 

yq  *  1 . 5 ,  ji  =  1 .  (a  be  ing  the  inner  wall  radius  and  b  being 

the  outer  wall  radius).  A  solid  curve  indicates  data  obtained  from 

a  numerical  solution  of  Eq.  (30);  the  dashed  line  is  a  plot  of  Eq. 

_3 

(31).  The  upper  pair  of  curves  is  for  v/yq  =  7.88  x  10  , 

2  2  -3 

w  /w  =  1.42,  the  lower  pair  is  for  v/ y  =  1.57  x  10  , 

p  o  o 

to  2/to  2  =  0.28. 

P  o 
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Fig.  15  Normalized  growth  rate  for  the  longitudinal  modes  versus  \>/ yq 

for  a  test  case  :  a=0.5m,  b=2.2m,  r^=0.99m,  i^l.Olm,  yQ=3,  £  =1- 
A  solid  curve  denotes  data  obtained  from  a  numerical  solution  of 
the  eigenvalue  problem  (30);  a  dashed  curve  denotes  data  from  the 
dispersion  relation,  Eq .  (31). 
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